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Abstract 

We report the performance of a 10 atm Xenon/trimethylamine time projection chamber (TPC) for the detection of X-rays (30 
keV) and 7 -rays (0.511-1.275 MeV) in conjunction with the accurate tracking of the associated electrons. When operated at such 
a high pressure and in 1%-admixtures, trimethylamine (TMA) endows Xenon with an extremely low electron diffusion (1.3 ± 0.13 
mm-cr (longitudinal), 0.8 ±0.15 mm-cr (transverse) along 1 m drift) besides forming a convenient ‘Penning-Fluorescent’ mixture. 
The TPC, that houses 1.1 kg of gas in its fiducial volume, operated continuously for 100 live-days in charge amplification mode. 



The readout was performed through the recently introduced microbulk Micromegas technology and the AFTER chip, providing a 
3D voxelization of 8mmx8mmx 1.2mm for approximately 10 cm/MeV-long electron tracks. Energy resolutions at full width half 
maximum (FWHM) inside the fiducial volume ranged from 14.6%(30keV) to 4.6%(1.275MeV). 

This work was developed as part of the R&D program of the NEXT collaboration for future detector upgrades in the search of the 
neutrino-less double beta decay (/3(30v) in 136 Xe, specifically those based on novel gas mixtures. Therefore we ultimately focus on 
the calorimetric and topological properties of the reconstructed MeV-electron tracks. In particular, the obtained energy resolution 
has been decomposed in its various contributions and improvements towards achieving the 1.4%^/lMeV/e levels obtained in small 
sensors are discussed. 


Key words: Double-beta decay, gamma and electron detection, microbulk Micromegas, time projection chamber, trimethylamine, high 
pressure Xenon, Fluorescent mixtures, Penning effect 
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1. Introduction 

The ability to perform 7-ray spectroscopy in conjunction 
with the accurate electron(s) reconstruction is of contempo¬ 
rary interest in fundamental science and technology. Such 
an asset is of utmost importance for instance in some imple¬ 
mentations of Compton cameras either in standard [1] or 
electron-tracking mode [2], 7 and X-ray polarimetry [3,4] as 
well as nuclear physics [5,6]. Generally speaking, a ‘dream’ 
scenario for detection is that where full event containment 
can be achieved in a medium with sparse ionization; an op¬ 
timal (‘intrinsic’) energy resolution emerges from the sub¬ 
sequent sub-Poisson fluctuations of the released charge [7], 
while the sense of the momentum vector can be obtained 
through resolving distinct track features such as character¬ 
istic X-ray emission at production and Bragg peak [8]. Re¬ 
construction of the track’s direction requires multiple scat¬ 
tering to be minimized and fidelity to the primary ioniza¬ 
tion trail to be ensured by proper design choices. If such is 
the case, accurate information about the photon energy, di¬ 
rection and sense of its momentum vector and polarization 
can be retrieved from the emerging lepton tracks, involv¬ 
ing in general multi-site cluster reconstruction algorithms. 
For very rare nuclear /3-decays, on the other hand, an ex¬ 
tremely high selectivity and background suppression can 
be achieved from the combined calorimetric and topologi¬ 
cal information [5]. 

2. High Pressure Xenon TPCs for MeV-y’s and 
electrons, and the case of NEXT 

High pressure Xenon time projection chambers (TPCs) 
capable of simultaneously reconstructing the energy and 
trajectory of MeV-electrons to great accuracy have been 
extensively discussed recently, and the reader is referred to 
[9,10] for details. As a gas, Xenon allows the initial ioniza¬ 
tion to span over a considerable distance even at high pres¬ 
sure (~10 cm/MeV at 10 bar) expanding the track details 
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and facilitating its reconstruction. At the same time, lep¬ 
ton tracks can be fully contained up to energies of the order 
of 10 MeV for a 100 kg/10 bar detector, by virtue of the 
high stopping power of Xenon. The scale of a cylindrical 
1 ton/10 bar detector with a height/diameter aspect ratio 
h/(j> = 1 approaches the asymptotic mean free path for pair 
production h = A e + e - = 9/7X 0 ~ 2 m, ensuring a high 7 
detection efficiency for all energies and event containment 
up to 10’s of MeV. 1 

The excellent tracking prospects of gaseous Xenon come 
at a modest x 2.5 worsening of the energy resolution relative 
to that of state of the art Germanium detectors: 


= 2.35 



0.45% 


IMeV 


( 1 ) 


I.e., for a fixed energy deposit (e), resolution is dictated 
by the Fano factor F (0.15 ± 0.02) and the average energy 
needed to create an electron-ion pair (Wj = 24.8 eV), [9]. 
In practical detectors this value can be approached through 
secondary light amplification (electroluminescence) [11], a 
process that benefits from pressurization [12]. For the si¬ 
multaneous reconstruction of the MeV-electron track, the 
separate optimization of the tracking (eg. with SiPMs) 
and calorimetric function (PMTs) has been pursued and 
demonstrated by the NEXT collaboration [13-15]. 

NEXT detection technique, based on the accurate re¬ 
construction of the 2-electron tracks emerging from the 
136 Xe /3/30i' decay with a total energy equaling Qpp = 
2.457 MeV, aims at suppressing the bulk of 7-ray back¬ 
grounds from natural radioactivity with rejection factors 
of ~ 10 7 , as required for the exploration of the inverted 
hierarchy of neutrino masses. At such a sensitivity level, 
achieving background-free conditions (if possible at all) re¬ 
quires advances on the topological analysis beyond the dou¬ 
ble Bragg-peak (double ‘blob’) identification long champi¬ 
oned by the Gottharcl collaboration [8]. Arguably, the de¬ 
tailed tracking in an axial magnetic field may provide a new 
path [16]: in this case 7-backgrounds display 1 blob and 1 
definite curvature (or chirality), while a /?/ 30is signal will 
double these features. For a relativistic track propagating 
over a distance l, the kick angle and the angular straggling 


1 X a is the radiation length. 




projected into the bending plane can be approximated, re¬ 
spectively, by [17]: 
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For practical magnetic-fields ( B) and pressures (P) this sets 
the desirable tracking step at the mm-level, below which 
multiple-scattering overrules magnetic bending. 2 Extract¬ 
ing additional topological information in the /3/30zz quest 
therefore implies that any blurring introduced by the dif¬ 
fusion of the initial ionization or the instrumental response 
should be minimized accordingly. Since this is not viable in 
pure Xenon, with an electron diffusion at the scale of 1 cm 
for 1 m drift [11], the possibility of a new generation of Xe- 
TPCs with a 10-fold increased track sharpness is being ex¬ 
plored. Achieving this goal, while preserving both the near- 
intrinsic calorimetric performance and the primary scintil¬ 
lation at the m-scale, presents a formidable challenge. 


3. Trimethylamine and its ‘magic’ mixtures with 
Xenon 


It has been noted that Xenon, when mixed with 
Trimethylamine (TMA: N(CH 3 ) 3 ) may, besides reducing 
the electron diffusion, form a mixture capable of simultane¬ 
ously displaying Penning effect and scintillation [18]. This 
concept is dubbed, hereafter, ‘Penning-Fluorescent’. The 
idea that Xenon VUV-light may be efficiently re-absorbed 
in fluorescent and ionization degrees of freedom of the addi¬ 
tive has appealing implications, allowing an improvement 
of the energy resolution through the reduction of the Fano 
factor and the increase of the ionization response 1 /Wj. 
as well as wavelength-shifting the primary and secondary 
scintillation. The possibility of using TMA for fine-tuning 
recombination in the search for directional dark-matter 
signals [19] as well as for charge neutralization of Ba ++ to 
Ba + (central to some Ba-tagging schemes, [16,20]) adds 
further value to this molecule, however these sophisticated 
scenarios (although qualitatively verified) remain specula¬ 
tive at the moment and are not discussed here. 

Trimethylamine is the smallest tertiary amine, and it is 
known to display a strong (and relatively fast) fluorescence 
from its lowest Rydberg 3s state (r 3s = 44 ns), with quan¬ 
tum yields up to 100% depending on the excitation wave¬ 
length \ exc . As for the case of ammonia, a large Stokes shift 
arises due to the transition between a pyramidal and planar 
geometry during the photon absorption-emission process. 
Hence, self-transparency can be anticipated for 1% admix¬ 
tures with noble gases up to 10’s of meters at 10 bar (Fig. 
1). Contrary to ammonia, the photo-absorption coefficient 


of TMA does not present marked vibrational bands, a con¬ 
sequence of the smearing introduced by the methyl rotor 
groups. This fact provides, on the other hand, additional 
electron cooling and so electron diffusion coefficients as low 
as Dy L = 250-350 /im/yTvm x Vbar have been reported 
in [21]. Such a diffusion would imply for a 10 bar detector 
a spread of the initial ionization at the scale of 1 mm-cr 
over aim path, as confirmed by this work, representing 
the smallest electron diffusion of any existing TPC that is 
known to the authors. 

The TMA fluorescence spectrum is centered at around 
300 nm and conveniently ends at the work function of cop¬ 
per (Fig. 1-down) thus reducing feedback from the photo- 
effect at metal surfaces. TMA itself is strongly opaque to 
Xenon light, with a mean free path of 0.4 mm for a 1% TMA 
admixture at 10 bar. 3 With additional VUV-quenching 
provided during the Xe-excimer formation stage (e.g., 
through the scintillating, Penning and photo-dissociative 
pathways made available by TMA) stable charge amplifi¬ 
cation up to several 1000’s can be achieved for pressures 
not exceeding 5 bar [21-23]. The gain enhancement rela¬ 
tive to that in pure Xenon has been often interpreted as a 
consequence of the close proximity of the TMA ionization 
potential to the energy of the excited states of Xe species, 
i.e., to the Penning effect. 
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Fig. 1. Top panel: World data compilation of TMA light properties at 
low pressure (1-100 mbar) for the photo-absorption coefficient, after 
Grosjean[24], Tannembaum[25] and Halpern[26]. Fluorescent quan¬ 
tum efficiency (QE) included as provided by Obi[27] (crosses) and 
Cureton[28] (circles), photo-ionization probability from Grosjean[24] 
(dot-dashed). Bottom panel: TMA and Xe scintillation spectra for 
different excitation mechanisms. The spectrum for 250 nm excitation 
[29] (QE=100%) is shown with thick line together with that for 193 
nm excitation as obtained by Cardoza in [30] (thin). Xenon scintilla¬ 
tion as measured by Koehler [31] for MeV-electrons at high pressure 
is indicated by dashed lines. Inset with molecule in 3D from NIST. 


2 ft is later shown that the spatial extent of an MeV-electron blob 
at 10 bar is 2 mm (1-cr). The suitability of mm-tracking will be 
demonstrated elsewhere, with the support of detailed simulations. 


3 Importantly, Xenon emission sits below the vertical ionization 
potential of TMA (Fig. 1), suppressing the track blurring that would 
be originated by photo-effect in the gas. 

























In this work it is described, mainly throughout section 4, 
the technical behaviour of a Xe-TMA TPC based on a spe¬ 
cific charge-amplifying structure (Micromegas). We later 
extend the experimental survey on electron swarm proper¬ 
ties of Xe-TMA mixtures at 10 bar in section 5, illustrating 
their impact in the reconstruction of extended tracks in sec¬ 
tion 6 . Section 7 is devoted to a discussion on the calorimet¬ 
ric response and the Penning-Fluorescent characteristics. 

4. Experimental setup: NEXT-MM 

4.1. System description and Micromegas readout plane 

NEXT-MM is a 73 liter cylindrical vessel made of stain¬ 
less steel, equipped with a field cage and 8 minx 8 mm pix- 
elated Micromegas anode, assembled and read out in TPC 
fashion. It has been developed by the NEXT collaboration 
to assess the performance of this family of micro-patterned 
charge readouts on large areas (700 cm 2 ) and high pres¬ 
sure, so as to accurately study MeV-electron tracks in low- 
diffusion mixtures for realistic conditions of event contain¬ 
ment. The TPC constituent materials are chosen to be ra- 
diopure: copper, PEEK, Teflon™ and polyimide, the main 
exception being the (solder-less) multi-pin signal connec¬ 
tors based on LCP . 4 The Micromegas were specially manu¬ 
factured at CERN with the ‘microbulk’ technique [32] pro¬ 
vided it grants a radioactivity content below 0.4 Bq/m 2 
for the 238 U and 232 Th natural chains [33]. Their progeny 
( 214 Bi and 208 T1) constitutes the main background in 136 Xe 
/3/30zx experiments [13,34]. Besides the presence of struc¬ 
tural copper, the TPC anode consists of a modest mass of 
3 g/m 2 of polyimide, from which much lower activities can 
be expected upon optimization of the fabrication process. 

A cross-section of the vessel together with a close-up of 
the anode region are shown in Fig. 2. The TPC has been 
designed to house about 1 kg of Xenon in its active volume 
when operated at 10 bar: h x </> = 38 cmx28 cm. A detailed 
description of the main TPC components, pressure, vacuum 
and recovery gas lines can be found in [22] and [35]. 

Technically, the TPC readout plane was made through 
a 4-fold sectorial arrangement with each sector enjoying 
the largest Micromegas manufacturable in the microbulk to 
date (~ 700/4 cm 2 ). The construction procedure relies on 
the chemical etching of doubly copper-clad polyimide foils, 
being performed in the present case on a 100 /jrn square pat¬ 
tern with 40 pm diameter holes placed at the square’s ver¬ 
texes. Within 50 /im of the holes’ entrance, corresponding 
to the thickness of the polyimide layer through which the 
gas amplification takes place, the second copper plate (later- 
acting as anode) is etched with a pattern of 8 mm x 8 nnn in¬ 
duction ‘pixels’. The number of amplification holes nears 7 
million. 


4 Liquid crystal polymer. ‘All-in-one’ sensor+cable designs combined 
with shielding is the most readily available solution to this persistent 
problem. 


For optimization of the Micromegas operation under 
Xenon-TMA we relied on ground-breaking work performed 
on small readouts in [22]. The best energy resolution there 
obtained was 1.4%-^/lMeV/s for 22 keV X-rays (FWHM), 
with a maximum operating gain of 400 at 10 bar. 



Fig. 2. Top panel: cross-section of the NEXT-MM vessel along the di¬ 
rection of gravity (downgoing) with its main elements being labelled. 
(T) Gas inlet; (2) gas outlet; (A) high voltage (HV) Teflon/copper 
feed-through; (4) Teflon cage and 241 Am source; (fT) PEEK bars for 
field-cage support; (6) insulating Cirlex™ rolled sheet; (7) pressure 
gauge; (5) Micromegas readout plane; ( 9 ) flat cable feed-throughs; 
@ TPC cathode; (lj) field cage rings; (12) vessel; (13) Nal(Tl) crys¬ 
tal and PMT; (M) 22 Na source. Bottom panel: up-view of the anode 
plane from about mid-drift. Uq) Footprint for signal out-routing; (10) 
after flat cable connection; (17) after connection reinforcement; (18) 
multi-pin SAMTEC connector before assembly; (19) Micromegas me¬ 
chanical reinforcement. The height of the field cage is h = 38 cm and 
the active diameter of the readout plane left by the rings’ shadow is 
<j> = 28 cm. 

4.2. Data taking and analysis 

Electron tracks used throughout this work were gener¬ 
ated by means of 7 -sources. The released primary ioniza¬ 
tion drifts towards the TPC anode where it avalanches in- 



side the Micromegas holes, inducing signals on its pixelated 
anode that can be recorded with specific sampling-ADC 
electronics, presently based on the AFTER chip [36]. In the 
absence of light sensors inside the TPC that would be capa¬ 
ble of recording the initial scintillation, the T 0 -information 
was provided by radiation emitted synchronously from the 
source, triggering the data acquisition system. For pres¬ 
sures around atmospheric the trigger signal was obtained 
with the help of a 241 Am source, the a-particles from its 
decays being conveniently tagged by a silicon diode in close 
proximity. Such a trigger system was encapsulated in a 
Teflon enclosure above the cathode and served as a source 
of quasi-simultaneous 237 Np X-rays in the range 25-60 keV, 
emitted through an opening in the cathode itself (Fig. 2- 
@). A 22 Na source located at about mid-drift was used 
for the measurements at 10 bar (Fig. 2-@, 7 -ray energies: 
0.511, 1.275 MeV). 

The experimental procedure started with a bake-out cy¬ 
cle of the vessel, gas injection and purification through re¬ 
circulation, and electric field conditioning. Typically, this 
lasted about 2 weeks. Following that and prior to the start 
of long data takings, the operating gain and drift field 
were systematically scanned for optimal working condi¬ 
tions: Edrift = 145 V/cm/bar, gain ~ 2000 (at 1 bar); 
Edrift = 80 V/cm/bar, gain ~ 200 (at 10 bar). Those were 
found to provide the best balance between Micromegas and 
field cage stability, highest signal-to-noise (S/N) ratio, con¬ 
venient drift times, low recombination and negligible at¬ 
tachment. It had been found in earlier works that TMA 
admixtures around 1% are convenient both in terms of Mi¬ 
cromegas gain response [ 22 ] and reduced diffusion [ 21 ] and 
so we concentrated on such a range. On the other hand, 
such concentrations would plausibly allow extracting the 
maximum light yield from the primary TMA scintillation 
as well as maximizing Penning (section 8 ). 

The data analysis proceeds starting from the 511-point 
waveform recorded by each of the 288 x 4 readout channels, 
with an event being a collection of waveforms. The analysis 
chain resembles, although supersedes, the one in [ 21 ]: 

(i) Zero-suppression. A global threshold is chosen rela¬ 
tive to the system equivalent noise charge (ENC) at 
about Eth = 5-10 xENC, set to minimize the presence 
of noisy channels and cross-talk: a) the non-firing pix¬ 
els and b) the non-firing time bins of the fired pixels 
are removed from the analysis chain. The threshold 
is defined over a certain charge pedestal, that is de¬ 
termined pulse-by-pulse and subtracted. 

(ii) Basic pulse shape analysis. Basic parameters of each 
waveform being above threshold are obtained: x-y 
pixel coordinates, position along the drift region (2), 
waveform integral (e), width ( <j t ), etc. 

(iii) Cluster finder. Two dimensional ionization clusters 
(projected in the x-y plane) are identified and tagged 
with a recursive depth-first search (DFS) method. 
For operational purposes a cluster is considered to be 
connected if a path exists between any two pixels be¬ 
longing to it; a path is a superposition of segments; 


keV 





Fig. 3. Several events obtained after calibration: x-y projection (left) 
and 2 -projection (right). Top panel: 1.275 MeV photo-electron. Mid¬ 
dle panel: 0.511 MeV photo-electron and characteristic X-ray emis¬ 
sion. Bottom panel: 0.511 MeV multi-Compton event. On the right 
column the thin lines represent the individual pulses registered for 
each pixel, and the thick line represents the event pulse (sum of all 
pulses). 

a segment is a line joining two neighbours; each pixel 
has 8 neighbours. Examples of events consisting of 1, 
2, and 3 clusters are given in Fig. 3. 

(iv) Identification of displaced X-rays. Displaced clusters 
with an energy in the range 20-40 keV are identified 
and undergo a Gaussian fit to refine their main pa¬ 
rameters relative to the estimate in (ii). 

(v) Quality cut. Events with fluctuating base-lines, 
anomalous noise or saturation are discarded. 

(vi) Volume fiducialization. Events not fully contained 
within 1 cm of the boundaries of the active region 
[</>, 0 -h\ are removed. 

(vii) Suppression of random coincidences. Random coin¬ 
cidences are suppressed by resorting to the physical 
correlation between a waveform’s width a t and its po¬ 
sition along the drift coordinate z. This correlation, 
stemming from the longitudinal electron diffusion, is 
particularly strong for X-ray clusters due their nearly 
point-like nature (see Fig. 10 for a graphical explana¬ 
tion). 

(viii) Event selection and analysis. Physical criteria are im¬ 
posed in order to select the events of interest. 

In the following, a magnitude A may be attributed to 
a generic event (A evt ), to a gamma-ray event (A 7 ), to an 
electron track within an event (A e ) or to an X-ray within 




















an event (A x ): A= x,y refers to energy-weighted positions 
projected into the readout plane, A = z refers to the arrival 
time of the maximum of the associated sum-pulse, A = a 
to its width, A = e to its sum-energy. 
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Fig. 4. Top panel: event position (energy-weighted) on the readout 
plane for a typical 22 Na run at 10 bar (with source placed at about 
the lower right corner, 5 cm distance from the readout boundaries). 
Micromegas readout sectors appear numbered according to their ge¬ 
ometrical position. Bottom panel: z event position (at the wave¬ 
form’s maximum) of a sub-sample of events displaying characteristic 
X-ray emission (thin line). The reconstructed position of the X-rays, 
smeared due to their finite mean free path, is shown with the thick 
line. The physical boundaries of the readout plane and drift region 
are indicated with continuous lines, while dashed lines indicate the 
fiducial volume. 


4.3. Calibration 


After these preliminary calibrations, a pixel-by-pixel gain 
equalization is performed iteratively by aligning the pixel 
energy distributions obtained over the x-y plane. Such pixel 
distributions are formed by histogramming the energy of X- 
ray clusters for which the contribution of the given pixel is 
the highest. This calibration procedure, that involves a mild 
sector-wise correction of gain transients (5-10%), converges 
after 2-5 iterations when the energy resolution reaches an 
optimum. 
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Fig. 5. Top panel: histogram showing the relative pixel gain (weighted 
by the number of events per pixel) before and after calibration, for 
1 bar (left) and 10 bar (right). Bottom panel: system resolution 
(FWHM) after calibration, suppression of random coincidences and 
TPC fiducialization for 1 bar (left) and 10 bar (right). Only the 
K a ,Kp lines and escape peaks are considered in the fit, with nor¬ 
malization free. The spectrum to the left is dominated by the escape 
peaks from the 241 Am 59.5 keV emission, while right shows charac¬ 
teristic Xe emission from the 22 Na 0.511, 1.275 MeV interactions. 


Data are first examined after step (v), Fig. 4. Reference 
positions are identified in the distribution of arrival times 
of the events’ ionization, providing an estimate of the z po¬ 
sition (Fig. 4-down). The slope of the cr t vs z correlation 
for X-ray clusters is stored for future cuts. Gain variations 
from sector to sector are corrected (10-20%) and the inte¬ 
gral of the waveforms converted to energy. 


The distribution of the per pixel ratios of the ~ 30 keV 
X-ray peak energy relative to the average value over the 
whole readout plane are shown in Fig. 5-up, each entry be¬ 
ing weighted by the number of X-rays counted in that pixel. 
These distributions, obtained before and after the calibra¬ 
tion procedure, inform about the robustness and quality of 
the algorithm, their final widths being below 1%(1 bar) and 


























2%(10 bar) FWHM, nowhere near influencing the system 
resolution obtained after calibration: 8.5±0.1%(1 bar) and 
14.6±0.2%(10 bar), Fig. 5-down. No sizeable deterioration 
of the latter values could be appreciated upon reduction of 
the statistics available for calibration by a factor x4. 

Fig. 6 shows the energy distributions and fitting func¬ 
tions used for calibration in a 4pixel x 4pixel region (squared 
area in Fig. 4), after calibration converged. A description 
of the fit details has been captioned. Although a peak- 
alignment with the anticipated ~ 1 % precision can be ap¬ 
preciated, there is a significant variability in the width of 
the main peak, beyond the statistical uncertainty of the fit. 




Fig. 6. X-ray energy distributions after calibration, obtained on a 
selected 4pixel X 4pixel region. Spectra are described through an it¬ 
erative multi-Gaussian fit consisting of K a , Kp characteristic Xenon 
emission and escape peaks. Peak positions and widths are bound rel¬ 
ative to the K a position by the known energy differences and 1/y/e 
scaling, with normalization left free. In case of bad fit quality or low 
statistics, a simple iterative gaussian fit is performed. Top panel: 1 
bar. Bottom panel: 10 bar. At high pressure only characteristic K a , 
Kp emission is expected, therefore the mean, width and normaliza¬ 
tion of the (unphysical) leftmost Gaussian are left free. This roughly 
accounts for the charge losses (due to a relatively high threshold 
around 5 keV) and increases the stability of the fit. 


4.4. Long-term stability and technical performance 

At 10 bar, about 11 million q-events from a 15 kBq 
22 Na source were taken during April-July 2014, totalling 
100 live-days. System operation at nominal working con¬ 
ditions proceeded continuously for 24 h a day and largely 
without manual intervention: no dedicated shifts were nec¬ 
essary except for sporadic monitoring of pressure and the 
current readings of the HV power supply for safety reasons. 
The bias of the readout plane was controlled with a N1471 
CAEN power supply, that would halt operation sector-wise 
in the event of an over-current exceeding 200 nA for more 
than 10 s. Effective data takings took place for about 30% 
of the total live-time, with the interval between acquisition 
runs dedicated to analysis and diagnostics. 

The level of pixel functionality was 90% throughout the 
measurements, the non-functionality approximately shared 
between permanently damaged pixels and connectivity is¬ 
sues (the latter discussed elsewhere [35]). Most of the dam¬ 
age arose during earlier measurements with a-sources and 
low-energy X-rays in [21,35] as a consequence of sustained 
over-current leading to violent sparking, the main causes 
identified being: i) readout mishandling, ii) a-emitter con¬ 
tamination of the cathode after prolonged use of a Rn 
source ([35]), and iii) damage during normal operation with 
7 -sources at relatively high gains. In case of damage, a 
manual intervention rendered the pixel electrically float¬ 
ing, restoring the electrostatic equilibrium in the affected 
region. The fringe fields resulting from this indirect bias¬ 
ing of the pixel potential do not penetrate beyond 1 % of 
the neighboring pixels, thus ensuring normal detector op¬ 
eration except for the modest geometrical loss associated 
to the intervention. 



time [days] 


Fig. 7. Evolution of the dark current and relative fraction of func¬ 
tional pixels in the Micromegas readout for 100 live-days under con¬ 
tinuous exposure of the TPC to a 15 kBq 22 Na source. The normal¬ 
ization is arbitrarily set to 1 at time=0, factoring out the damage 
occurring during commissioning (5%) and the presence of uncon¬ 
nected pixels (see text for details). Three pixels were damaged dur¬ 
ing this prolonged operation (sector 2 and 4). Grey rectangles show 
the periods when data takings took place. 

Although dark currents normally stayed within 0.015- 
0.03 nA/cm 2 , values as high as 0.3 and even 3 nA/cm 2 were 
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reached in two of the sectors (Fig. 7). This did not affect 
detector operation except for a somewhat x 1/2 lower gain 
in the most extreme case. The absence of a strong gain 
drop or hot spot suggests that these leaks may be related to 
loss of insulation downstream the amplification gap. Weak 
points may appear at connectors and cable traces that re¬ 
sult indirectly biased when a pixel is left floating during the 
aforementioned process. 

The presence of micro-discharges in non-insulated Mi- 
cromegas operated close to breakdown is a well documented 
phenomenon (see for instance [37] and references therein), 
resulting in the present case in an average rate of current 
excursions about 10 /min/m. Provided the pre-amplifier 
input is protected they have no influence in the detector 
stability. As shown below, they have also a negligible im¬ 
pact in the system exposure, since their natural scale is that 
of few amplification holes. In most cases this phenomenon 
affects only one pixel (C p i x = 20 pF) during a short period 
of time. More intense current excursions are exponentially 
less frequent although potentially more harmful. 

Four types of discharges can be found in this system, 
qualitatively enumerated below in order of increasing 
strength: 

I Low intensity: current excursions at the level of 20-60 
nA. 

II Medium intensity: current excursions exceeding the 
HV current-limit of 200 nA, resulting in a HV trip. 

III High intensity: current excursions exceeding the HV 
current-limit of 200 nA during 10 s, resulting in a HV 
ramp down. 

IV Permanent damage: a high intensity discharge that 
permanently develops a high conductivity path. 

Incidentally, stronger micro-discharges were found to 
have also a stronger impact on the fraction of system ex¬ 
posure that is affected, This magnitude can be 

approximated by the micro-discharge rate, and the area 
( A a ff ) and time (A T a ff) affected after each occurrence, 
as detailed in Table 1: 

~ rateAA aff AT aff (4) 

After accounting for all instabilities in the readout plane, 
the overall reduction in exposure is at the level of 0.5 x 10~ 3 , 
seemingly of little practical importance. Moreover, if an au¬ 
tomatic method would be implemented for leaving a dam¬ 
aged pixel in open-circuit after IV-type discharges, opera¬ 
tion without manual intervention becomes possible in nr 2 
experiments for Micromegas operating close to breakdown. 

A priori more worrisome for long-term operation is the 
initial damage and the rate of new damage. As for the ini¬ 
tial 5% damage, it is difficult to state which part is un¬ 
avoidable during the commissioning and system condition¬ 
ing. Certainly 1% levels would be approachable through a 
factor x 2 reduction in pixel size and improved system han¬ 
dling. Implementation of a tight quality control (QC) would 
bring additional gain. The chosen TPC orientation maxi¬ 
mizes the accumulation of dust, that should be avoided in 


type 

rate 

AT V/ 
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A A-Ioss 
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I 

10 /(min m 2 ) 

2 s 

pixel 

- 

0.02 x 10 -3 

II 

2 /(h m 2 ) 

10 s 

sector 

- 

0.09 x 10 -3 

III 

5/(d m 2 ) 

2 min 

sector 

- 

0.14 x 10“ 3 

IV 

i%/y 

2 h 

sector 

pixel 

CO 

1 

o 

X 

o 

CQ 
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Table 1 

Types of micro-discharges in the Micromegas readout together with 
the (approximate) affected time, affected area and damage, if any. 
In case IV the time is merely indicative since it involves manual 
intervention. The last column shows the fraction of experiment’s 
exposure affected, amounting in total to less than 0.05%. 

an experiment if no particular geometrical constraint ex¬ 
ists. The remaining damage rate of 1%/year would benefit 
in a similar way from the above measures. 

5. Response to 30 keV X-rays and comparison 
with earlier results 

5.1. Energy resolution and intrinsic performance 

The energy resolution obtained for X-rays after calibra¬ 
tion is compiled as a function of the drift field in Fig. 8 (cir¬ 
cles). The horizontal bands show, as a reference, the energy 
resolution obtained in small setups and for similar gains 
[22]. We refer to this resolution, hereafter, as ‘intrinsic’ (to 
the charge-amplification readout). Since such an intrinsic 
response offers a promising extrapolation to the Qpp en¬ 
ergy we discuss below the observed excess. 
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Fig. 8. Measured energy resolution (FWHM) as a function of the drift 
field in the TPC fiducial region (circles), together with trend lines. 
Noise-subtracted values are given by the dashed line, while squares 
indicate the result after both pixel selection and noise subtraction, 
for the highest fields. Shadowed bands represent the values obtained 
earlier in small setups and for similar gains in [22]. 

On the one hand the degradation with lower drift field 
cannot be attributed to attachment or recombination since 
charge losses from 80 to 40 V/cm/bar amount to some 10% 
(subsection 5.2), hence implying a negligible contribution 
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to the resolution below 3% (to be summed in quadrature). 5 
Instead, the reason for the observed behaviour can be found 
on a relatively low S/N ratio and a reduced sampling fre¬ 
quency at low drift fields. The latter is determined by the 
limited number of sampling points provided by the elec¬ 
tronics (511) in connection with the reduced drift velocity. 
In fact, given the response function of the AFTER chip for 
5-excitation: 

f(t) = Ae~ 3t/r (t/r) 3 sin(t/r) (5) 

even for the longest shaping time available (r = 2/zs), the 
minimum sampling time achievable (if recording the com¬ 
plete drift region) barely satisfies the Shannon-Nyquist con¬ 
dition for the lowest drift fields. Additional contributions 
to this response such as the ion transit time along the Mi- 
cromegas amplification gap or the extension of the electron 
track are expected to be small and indeed r = 2.1/is de¬ 
scribed well the temporal width of X-rays found close to 
the anode. 

From the above perspective it can be expected that the 
deterioration of the energy resolution at low fields comes 
from the inability to accurately reconstruct the integral of 
f(t), as illustrated by simulation in Fig. 9-up. For clarity, 
f{t) displays here the amplitude of a 30 keV charge deposit 
after including the longitudinal charge spread at about mid¬ 
drift (black line), as seen by the acquisition system. The 
resulting waveform is randomly sampled with the corre¬ 
sponding time bin A Tbi n and smeared with white noise, 
a Gaussian fit finally implemented exactly as done in the 
data analysis (grey lines). Simulations indicate that, if as¬ 
suming ENC= 7500e~ (50% in excess of the average ENC 
measured over the whole readout plane), most of the ob¬ 
served deterioration in the reconstructed pulse integral can 
be explained (dashed lines in Fig. 8). 

A second source of deterioration can be found in the 
aforementioned variation of the energy resolution from 
pixel to pixel, illustrated in Fig. 9-down by histogramming 
the energy resolution of those pixels with higher statistics. 
Even if accepting the presence of statistical and systematic 
effects in the fits of Fig. 6, some additional contribution 
seems to be needed to explain the spread of these distri¬ 
butions. In particular, the larger spread observed at high 
pressure can be qualitatively understood from the steeper 
behaviour of the gain vs amplification field (E a ) systemat- 
ics in [22], with the quantity A g increasing with 
operating pressure. A recent microscopic simulation of the 
Micromegas readout shows that the energy resolution at 
high pressure is indeed much more vulnerable to the ge¬ 
ometrical tolerances in the diameter of the amplification 
holes, by a factor x3 [41]. If approximately selecting the 
best 10% pixels (outliers removed) values closer to intrinsic 
can be found (shadowed bands in Fig. 9). An approximate 


5 This follows the canonical formula for the case where recombi- 
nation 7 Z is decoupled from other processes: 2.35^/72./(1 — 7 Z) n 0 , 
n 0 = e/Wi ~ 1200e~ for 30 keV X-rays in Xenon. 



lbar: (events per pixel>200), 1 Obar: (events per pixel>500) 



Fig. 9. Contributions to the energy resolution beyond the intrinsic 
expectation from the Micromegas readout. Top panel: simulation of 
the reconstruction algorithm (grey lines) for the assumed response 
function given in black at the nominal thresholds (dashed). Left and 
right illustrate typical running conditions for low and high pressure, 
respectively. Bottom panel: distribution of the pixels’ energy resolu¬ 
tions (FWHM) for 1 and 10 bar. Shadowed bands correspond to an 
approximate selection of the best 10% after removing outliers. 


decomposition of the energy resolution for the highest drift 
fields studied is: 


9W= ^0.075 2 (int)+0.002 2 (S/N)+0.035 2 (p-p) (6) 

$Hi 0 fcar=v / 0-09 2 (int) + 0.07 2 (S/N)+0.09 2 (p-p) (7) 

that includes the contributions intrinsic to the readout 
(int), signal to noise (S/N) and pixel-to-pixel (p-p). It must 
be noted that collimation at the level of few mm 2 (much 
smaller than the pixel size in NEXT-MM) is often used in 
test setups and that there is an implicit selection bias in that 
under-performing wafers are discarded before characteri¬ 
zation or suffer from early damage. Microbulk Micromegas 
even on modest sizes of 10’s of cm 2 as those used in CAST 
already show energy resolutions mildly higher than intrin¬ 
sic [38]. Put simply, any residual gain variation below the 
minimum resolvable unit of a pixel (or a strip, in that case) 
is hardly correctable during an experiment. 







5.2. Parameters of the electron swarm 

The precision in the topological reconstruction of the ion¬ 
ization trail left by a primary track in a TPC depends on 
the native system voxelization (pixel size, time bin) and the 
‘swarm’ properties of the electrons released: drift velocity 
Vd and pressure-normalized diffusion coefficients D* L T . The 
reconstructed energy may be limited, on the other hand, 
by losses due to electron-ion recombination or electron at¬ 
tachment (1Z, 77): 


Qe{z) _ vz 

Qoo(0) 

(8) 

&z,xy(z) = D lt -j= 

(9) 

= — 

Vd 

(10) 


Qe(z) refers to the charge collected as a function of the 
electric field and drift distance, and Qoo(0) represents the 
amount of charge released by the event in the absence of re¬ 
combination. These electron swarm parameters can be ex¬ 
tracted from the same X-ray sub-sample used for the calori¬ 
metric analysis of subsection 5.1, provided their point-like 
nature allows approaching the d-type ionization conditions 
assumed in eqs. 8-10: i) the drift velocity is determined from 
the difference in arrival times of X-ray clusters produced at 
two well defined positions and ii) the longitudinal diffusion 
D* l and attachment coefficients are obtained from the of vs 
0 and e vs 2 cluster behaviour (eqs. 8, 10). It is in principle 
possible to extract the transverse diffusion coefficient Dif 
through a separate analysis, however the spatial sensitivity 
at 10 bar was found not to be sufficient for this. Details on 
the extraction of D* L T , Vd for 1-3 bar can be found in [22] 
and so here only the 10 bar analysis is described. 

The drift velocity was determined from the TPC anode 
and 22 Na source positions, the latter identified as a pro¬ 
nounced peak in the X-ray time distribution, correspond¬ 
ing to z = 19 ± 2 cm (Fig. 4). Arrival times from clusters 
in the TPC-anode region are synchronous with the Nal(Tl) 
trigger signal and could be identified with an estimated 2 
cm uncertainty. 

Concerning longitudinal diffusion, the of vs z correlation 
considerably widened at 10 bar as compared to 1 bar (Fig. 
10). This is consistent with the 10-fold S/N reduction and 
the associated deterioration in the reconstruction of the 
pulse characteristics. In order to minimize the additional 
blurring introduced by random coincidences at 10 bar, the 
main electron track associated to the X-ray was requested 
to be in the range z e = 19 ± 5 cm. Linear fits in a selected 
region in z around the source position yielded the desired 
slopes (shown in Fig. 11). For illustration, the simulation 
of the pulse reconstruction procedure at various z-positions 
was performed, and the final lineal fit is shown for the 80 
V/cm/bar case (dashed line). The parameters used in the 
simulation are r = 2.15 /rs, ENC=7500e“, A Tbi n = 1 /xs, 
Sth = lOxENC, gain= 200, and the slope introduced was 
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Fig. 10. Temporal variance of the pulse as a function of the drift 
position for isolated X-rays. The boundaries of the drift region as 
well as the trends are indicated by continuous lines. Dashed lines 
together with the drift boundaries delimit the cut region selected for 
suppressing random coincidences. Top panel: 1 bar (source placed 
above cathode). Bottom panel: 10 bar (source placed at mid-drift). 
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Fig. 11. Determination of the slope from the cr^ vs 2 : plots for X-ray 
clusters. An iterative Gaussian fit has been performed in 1 cm steps 
and the slope obtained from a weighted fit. The plot at 80 V/cm/bar 
includes the results from simulation (dashed lines). 

the measured one. It is important to note that the system- 















atic upwards shift (resulting from the effective increase of 
the pulse width during reconstruction) is correctly captured 
in simulation, and it is qualitatively reproduced for lower 
fields where the effect is considerably larger. However, even 
in this highly non-ideal situation, the reconstructed slope 
stays within 5% of the value introduced, at most. The un¬ 
certainty in the determination of the longitudinal diffusion 
coefficient is obtained from this 5%, the 10% uncertainty 
in the drift velocity, the statistical error of the linear fit 
and the variations for different assumptions on the fitting 
range. When summed in quadrature, the resulting uncer¬ 
tainty lies in the range 10-15%. 

Fig. 12 shows a compilation of the main results in the 
field region that is of interest for the NEXT experiment to¬ 
gether with a comparison with the microscopic code Mag- 
boltz 10.0.1 [39]. Additional measurements performed in 
small setups at various pressures can be found in [21] (drift 
velocity) and [19] (drift velocity and longitudinal diffusion 
coefficient), showing good mutual agreement. The com¬ 
bined results convey evidence of the strong electron cool¬ 
ing properties of TMA, showing a reduction in diffusion 
within factors of x3(L), x20(T) of the ones achievable in 
pure Xenon, and significantly outperforming the expecta¬ 
tion for Xe-CH 4 admixtures [8] by nearly a factor x2. Xe- 
TMA sets the diffusion scale in Xenon mixtures to 1 mm 
for 1 m drift at 10 bar, less than a part in 200 relative to 
the size of a Qpp event at the same pressure. 

Interestingly, the cross-over shown by the drift velocity 
in Xe-TMA mixtures relative to pure Xenon is well repro¬ 
duced, and so below 75 V/cm/bar the mixture becomes ac¬ 
tually slower than pure Xenon. Since this type of systems 
often operate in a closed gas loop, a high sensitivity of the 
drift velocity to the reduced field may turn to be a rele¬ 
vant practical advantage, allowing the optimization of the 
voxelization along z during operation. 

The dependence of the X-ray energy peak with the drift 
distance is shown in Fig. 13-up for Edrift = 80 V/cm/bar. 
Further to the z e cut, random coincidences have been sup¬ 
pressed through a of-z correlation cut (dashed lines in Fig. 
10). The resulting e vs 0 distribution is binned and fitted 
iteratively to a Gaussian function, yielding for the high¬ 
est field an attachment coefficient rj = 7% ± 4%/m after 
an exponential fit (corresponding to an electron lifetime of 
12 ms). 

At last, the relative fraction of primary ionization Q that 
is recombined can be obtained by studying the variations of 
the X-ray energy peak position (back-converted to charge), 
obtained at about mid-drift and corrected for attachment. 
Due to the anticipated smallness of the effect the analysis 
has several caveats, in particular in view of the variation of 
the sampling time and diffusion with field. Simulation in¬ 
dicates that any systematic shift of the reconstructed pulse 
integral is below 1% even for the lowest fields, and diffusion 
is relatively stable (Fig. 12), so we may neglect any effect re¬ 
lated to the pulse shape reconstruction. The remaining sys¬ 
tematic effects are connected to the presence of (not under¬ 
stood) gain transients at the level of 5% and to the stronger 
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Fig. 12. Compilation of the main properties of high pressure Xenon 
TPCs for 1(1.0) 3(2.7) and 10(10.1) bar and different TMA admix¬ 
tures, together with simulation (Magboltz 10.0.1, [39]), as a function 
of the reduced electric field. Left: drift velocities. Right: longitudinal 
and transverse diffusion coefficients. Simulation results for Xe-CH 4 
(Gotthard) and Xe-TMA in a 0.1% admixture are shown for illus¬ 
tration. 
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Fig. 13. Top panel: X-ray energy peak position as a function of drift 
distance for Edrift = 80 V/cm/bar. Bottom panel: X-ray energy 
peak position (back-converted to charge) relative to 80 V/cm/bar 
after correcting for attachment (only statistical errors shown). 

^-dependence of the reconstructed energy at low fields, im¬ 
plying up to a 10% correction. Under these assumptions a 
coarse estimate of the maximum recombination observed 
in present X-ray data is 1Z X = 1 — Q4o(0)/Qso(0) = 0.12 ± 
0 .02(sta)tg os(sys), the statistical uncertainty coming from 
sector-to-sector variations (Fig. 13-down). Within such a 
large uncertainty, this value is compatible with the charge 
loss reported in [41], that sets a scale for recombination at 
the level H x = 1 - Q 8 o(0)/Qoo(0) < 0.20 in 1% TMA ad¬ 
mixtures at 10 bar. 
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6 . Extended electron tracks from 7 -rays 

A typical 7 -ray coincidence spectrum obtained for single¬ 
tracks after calibration is shown in Fig. 14. The TPC 
working conditions are: Xenon/trimethylamine filling gas 
at 99/1 volume admixture, P = 10.1 bar (T = 20 deg), 
Edrift = 80 V/cnr/bar, E amp = 11.1 kV/cm/bar. These 
settings correspond to a gain for the Micromegas readout 
of about x200, with the main TPC charge-drift properties 
being Vd = 1.16 ± 0.12 mm//is, a L (lm)= 1.3 ± 0.13 mm, 
<r T (lm)= 0.8 ± 0.15 mm (from Fig. 12), 77 = 7% ± 4%/rn 
(Fig. 13). Indirect estimates for initial charge recombi¬ 
nation, Penning transfer probability and primary scin¬ 
tillation yield are 1Z X < 0.2 (section 5.2), rp > 0.1 and 
S'i(~300nm) > 50 ph./MeV (given later in section 7.2), 
respectively. The remaining parameters impacting the ac¬ 
curacy of the reconstructed track are the threshold energy 
£th = 4 keV (ENC=5000e _ ), the response function (dom¬ 
inated by the AFTER chip) of 1.2 /rs-cr (1.4 mm) and a 
time sampling of 1 /ts ( 1.2 mm), that ultimately results in 
a native 3D-voxelization of 8 nnnx 8 mmxl.2 mm. 

single track spectrum 



energy of the event [keV] 


Fig. 14. Single-track energy reconstruction in the TPC after cali¬ 
bration (thick line), for events coming from a 22 Na source in coin¬ 
cidence with a Na(Tl) detector. The contributions from tracks fully 
contained in each of the individual sectors ( 61 , 2 , 3 , 4 ) are also shown. 

The most prominent features in the spectrum correspond 
to the photo-peaks from the annihilation of the e + emitted 
in the 22 Na decay, at 0.511 MeV, and the 7 -emission after 
the daughter’s decay at 1.275 MeV. The energy resolution is 
not sufficient to separate the escape peak associated to the 
0.511 MeV interactions, appearing as a shoulder to the left 
of the main peak. For sector 2 (where the gain was lower) 
and sector 1 (that showed larger gain transients) the tails 


to the left are more pronounced and the resolution wors¬ 
ened. The suppression of the 0.511 MeV region in sector 
3 relative to the high energy Compton contribution is due 
to the trigger geometry (the solid angle for 0.511 MeV is 
determined by the coincidence with the Nal(Tl) detector, 
covering mostly sectors 1, 2 and 4). 

6.1. Calorimetry 

We focused on events that satisfy the analysis criteria v- 
vi (subsection 4.2), and concentrate in particular on elec¬ 
tron tracks displaying characteristic X-ray emission. This 
additional criterion allows suppressing the escape peaks 
and multi-Compton events as well as random coincidences. 
The latter can be achieved through the er/ vs z correlation 
cut on the X-ray cluster. 6 

An example of a typical 0.511 MeV event is shown in 
Fig. 3 (middle panel). Electrons with this energy deposit 
about 1/2 of their energy in the track’s last stage as a re¬ 
sult of the Bethe-Bloch increase (Bragg peak) and the in¬ 
creased multiple scattering, a region generically referred to 
as ‘blob’. The approximately equal energy sharing between 
the blob and the (nearly) ‘mip’ region results in relatively 
small tracks of 3-4 cm length. As later shown, electrons 
from the 1.275 MeV region can reach, on the other hand, 
considerably longer extensions of up to 10-15 cm (Fig. 20). 
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Fig. 15. Left: energy distribution around 0.511 MeV obtained for 
double-tracks (i.e., events consisting of a single-track in the presence 
of a displaced X-ray), shown for the TPC fiducial region and for sector 
4. Right: energy distribution around 1.275 MeV for single-tracks. 

The spectrum obtained in the 0.511 MeV region for the 
full TPC after fiducialization (1.1 kg) is consistent with 
an energy resolution of 8.5 ± 0.1% (Fig. 15-left), while the 
sector with the highest statistics and S/N (S 4 ) performs 
at 7.0 ± 0.1%. The energy resolution obtained in the 1.275 

6 It can be noted that, coincidentally, the selected sub-sample is 
the one used for extracting the electron swarm parameters, where 
emphasis was put on the X-ray behavior. 
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MeV region can be estimated from the high statistics 
single-track sample as 4.6±1.5%(TPC) and 3.9±1.0%(S4), 
Fig. 15-riglit. Therefore we come to energy resolutions for 
MeV-electron tracks of 4.5-5.5%y / lMeV/e (0.511 MeV) 
and 4.4-5.2%y / lMeV/e (1.275 MeV), compatible with 
l/y^-scaling. These numbers are within a factor x2 of 
the 2.5%y/lMeV/e obtained for X-rays in section 5.1, that 
deviates from the intrinsic energy resolution of the Mi- 
cronregas readout at 10 bar: 1.4%^/lMeV/e. A plausible 
explanation for the deterioration of the calorimetric re¬ 
sponse to X-rays has been attributed to the low S/N and 
non-uniformity of the readout. However, the energy reso¬ 
lution obtained for MeV-electron tracks suffers from yet 
another x2, for which it seems that a different source is 
needed. 
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Fig. 16. Position of the 0.511-MeV peak as a function of the drift 
distance for 80 V/cm/bar and 10 bar. 

It must be noted that the positions of the 0.511 MeV and 
30 keV energy peaks are well reproduced after X-ray cal¬ 
ibration within 2% (Figs. 13, 16), showing little evidence 
for attachment or strong additional recombination besides 
that observed for X-rays. A marginal correlation of the re¬ 
constructed energy with the track orientation and the num¬ 
ber of pixels fired was observed, that may be an indica¬ 
tion of an imperfect calibration, however it was just too 
small to be easily correctable, and clearly does not account 
for the observed deterioration. Very modest, < 0.5%, im¬ 
provements could be achieved by selecting specific regions 
along or across the chamber. The main explanation in hand 
comes from a relatively large system cross-talk originated 
during signal transmission (1/20-1/30), that is character¬ 
ized in the present case by perfectly unipolar signals, and 
that appears synchronized with the event’s largest pulses. 7 
A soft anti-coincidence condition in a ±4/zs (±5nnn) win¬ 
dow around the time of the event allows to suppress them 
when they appear isolated, however it is virtually impossi¬ 
ble to avoid their contribution if they are produced onto the 
pixels already having a signal above or just near threshold. 


7 The best explanation for this anomalous cross-talk is the presence 
of < 1 OkQ resistive paths or dielectric shunting conductances be¬ 
tween some of the traces (either at the cables or at the footprints), 
something that had been qualitatively identified earlier. The com¬ 
plexity associated to the unambiguous identification of a cross-talk 
event, due to the presence of net charge, discouraged us from the 
production and study of detailed cross-talk maps. 


A ball-park estimate treating cross-talk as correlated noise 
with a flat probability in the range 0-1/20 would yield 3.5% 
additional peak broadening, while an average of 1.7 pix¬ 
els showing cross-talk would plausibly explain the observed 
discrepancy between X-rays and extended electron tracks. 

6.2. Track topology 

Besides the calorimetric function, an essential feature for 
f3(3 applications is the accurate reconstruction of the elec¬ 
tron track end-point (blob). This is illustrated in Fig. 20 
through several hand-picked single-track events in the en¬ 
ergy region around 1.275 MeV. In this (biased-by-selection) 
sub-sample, that roughly describes the behaviour of about 
1/3 of the events in this energy region, the end-blob is vis¬ 
ible to the naked eye. The remaining 2/3 consist of cosmic 
rays, highly bent tracks and multi-Compton events. The 
selected energy approaches the most probable electron en¬ 
ergy in the neutrino-less decay mode (Qpp/2), hence effec¬ 
tively indicating how ‘half’ of the expected signal from the 
long sought decay looks like. 

The main practical limitation for recognizing a double¬ 
blob f3f30is event is the inability to adequately ‘un-fold’ 
the tortuous electron tracks into a straight line where the 
blob(s) cannot be mistaken. If no attempt is made to quan¬ 
tify the final figure of background/signal suppression, the 
problem can be analyzed through a sub-sample of tracks 
selected to be straight along a particular spatial direction. 
As will be immediately shown, 0.511 MeV tracks are long 
enough to allow a separation of the nrip and blob region, 
and so their analysis universally describes the final part of 
any electron of arbitrary energy in the present experimen¬ 
tal conditions. Provided the system’s native voxelization 
is 8rnmx8nnnx 1.2mm we focus on the z coordinate due 
to its higher resolution: Fig. 17 shows the superposition of 
all event pulses, aligned with respect to the position of its 
maximum and divided by the number of events, for differ¬ 
ent track lengths along the ^-coordinate. 

As expected, the Bragg peak can be distinctly identified 
in extended tracks, with a blob to nrip energy ratio of 5. 
The blob width is 2.5 mm-a, being only partly explained 
by the electronics response function (1.4 nun), and with 
negligible contributions from binning (0.3 nnn) and diffu¬ 
sion (0.55 nnn). This sets the physical scale of the electron 
blob at about 2 mm-cr (for 10 bar Xenon) implying that in 
present conditions not just the blob but its details start to 
be resolvable. 

The influence of the pixelization in track reconstruction 
can be easily appraised by performing the same analysis 
for the longest tracks in all 3 directions of space (Fig. 18). 
Clearly, the wider pixelization (8 mm) results in a loss of 
information, but how important is it?: Fig. 18 describes the 
average behaviour and extension of the final energy blob in 
MeV-electron tracks but does not give an idea about the 
energy loss fluctuations. By approximately selecting 8nnn- 
regions both around the ‘blob’ and ‘nrip’ regions (bands in 
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Fig. 17. Average energy loss profile for 0.511 MeV-electron tracks of 
different lengths along the 2 direction (1.2 mm bin). Continuous and 
dashed lines represent upgoing and downgoing electrons, respectively. 




Fig. 19. Event-by-event blob and mip energy distributions obtained 
from a selection of ‘long’ tracks (i.e. 3.2cm< lengthy z < 4.8cm), 
by selecting different Ax , y , 2 bands in each projection. Up: 
A x,y = 8mm (left) and A 2 = 8.4mm (right). In the left figure the 
spectra corresponding to tracks extended in x are indicated by contin¬ 
uous lines, whereas they are dashed for y. Bottom-left: A 2 = 3.6mm. 
Bottom-right: native 2 -binning. 

interpret a 8 mm mip region as a blob is below 1/20 at a 
signal efficiency of 96%, if a sharp cut at 150 keV is used. 
The situation improves for smaller Az intervals until the 
mip and blob distributions completely decouple if the na¬ 
tive z pixelization (1.2 nun) is chosen. This again suggests 
that the blob region is being resolved in present experimen¬ 
tal conditions and that nun-scale tracking capabilities are 
an important asset. 


Fig. 18. Average energy loss profile for 0.511 MeV tracks selected to 
be long (3.2-4.8cm) along any of the 3 spatial projections. Binning 
is 8 mm for x, y projections and 1.2 mm for 2 , following the native 
TPC voxelization. Continuous/dashed lines represent (from top to 
bottom) right/left, forward/backward, up/down -going tracks. Bands 
implement the definition of blob and mip region used in text. 


Fig. 18) identification can be attempted track-by-track, re¬ 
sulting in the distributions shown in Fig. 19-up for x and 
y (up-left) and z (up-riglit). The probability to mistakenly 
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Fig. 20. 3D reconstruction and 2D projections of several hand-picked electron tracks in a ±50 keV region around the 1.275 MeV peak, shown 
with a voxelization of 8 mm x 8 mm x 10 mm (x, y , z). Upper 2 figures for each event share color bar. 





























































































































strated. 


7. Discussion 

7.1. System aspects 

From the calorimetric point of view the performance of 
NEXT-MM lkg/lObar Xenon TPC compares favourably 
with earlier results for similar systems (Table 2) with ex¬ 
ception of the NEXT baseline technology, based on electro¬ 
luminescence. The energy resolution for a charge-readout 
system may be decomposed according to the contributions 
that have been identified in this work as: 

^ ^-^^aLt+vl-p + Pl/N+vl-talk+vlal ( U ) 

^ F+ f +n \j (1 WI n)e ( 12 ) 

From left to right the terms in eq. 11 correspond to the in¬ 
trinsic contribution from the Micromegas readout (includ¬ 
ing effects in drift and amplification regions), gain spread in 
the amplification holes (experimentally manifested as pixel- 
to-pixel variations of the resolution after gain equalization), 
limited S/N and sampling frequency, cross-talk, and cali¬ 
bration effects. Eq. 12 shows the canonical formula for the 
intrinsic energy resolution under the assumption of decou¬ 
pled fluctuations, with / corresponding to the variance to 
gain-squared ratio of the multiplication process. Once / 
is obtained from microscopic simulation, this formula de¬ 
scribes small Micromegas X-ray data to good accuracy [41]. 
Additional physical limitations to eq. 12 for MeV-electron 
tracks can be reasonably discarded: i) Recombination for 
0.511 MeV tracks is from present data compatible with that 
observed for X-rays 'R-^n 1 /lZ x ~ 0.98 ± 0.02 while a small 
effect is present at 1.275 MeV 7 ^ 12757 /^sivy = 0-97 ± 0.01; 
ii) physical characteristics of the track like bremsstrahlung 
[9] are determined by the amount of Xenon and therefore 
any contribution to the energy resolution is largely con¬ 
strained by the most accurate measurements in table 2; iii) 
effects related to the track orientation and charge density 
in the amplification holes were too small to be isolated in 
this analysis. 

From the previous arguments it may be assumed that the 
measured energy resolution is dominated by instrumental 
limitations and not by the different physical behaviour of 
X-ray and 7 -tracks. Several sources of improvement can be 
identified: i) according to [36] noise minimization down to 
ENC= 2000e~ is feasible with a reduced cable capacitance 
down to 50 pF {C P i X +C ca bi e — 180 pF in this system [35]) 
while increased sampling frequency (x 4) and cross-talk op¬ 
timization in cables and connections should suppress the 
most important contributions to eq. 11 : <Js/n and &x-talk- 
The implementation of well targeted QC procedures in con¬ 
nection with hole-pattern optimization, geometrical gain 
compensation [42] or increased fabrication accuracy could 
bring further improvements in the magnitude and unifor¬ 
mity of the operating gain. This will help at further reduc¬ 
ing the terms as/N and cx p _ p , but it remains to be demon- 
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Table 2 


Energy resolution (l/y^-extrapolated to 1 MeV for comparison) 
in Xenon-based detectors and for different electron energies. For 
reference, values measured in state of the art microbulk Micromegas 
under localized irradiation and the Fano factor limit are shown. 
They approximately represent intrinsic limits to charge-amplification 
readouts and electroluminescence, respectively. 

Besides the good energy resolution observed for 0.511 
MeV electrons these type of readouts have remarkable 
strengths, like the easiness and scalability of the readout 
segmentation, simplicity and high fidelity to primary ion¬ 
ization. The chosen voxelization of 8 mmx 8 nnnxl.2 mm 
provides a clean separation between blob and mip regions, 
largely unaffected by diffusion even at the lm-scale. The 
detector can work stable for 100 live-days without signs of 
deterioration except for a mild damage rate extrapolable 
to about 1 %/year or less due to conditioning effects. 

7.2. Penning-Fluorescent behaviour 

Indirect evidence for Penning-transfer reactions in Xe- 
TMA mixtures has been obtained in wire-chambers [23], 
GEMs [44], and Micromegas [22]. Their scintillating prop¬ 
erties on the other hand have been reported recently [40]. 
If considering only energy transfers between the lowest ly¬ 
ing states, the behaviour of a (VUV-quenclied) Penning- 
Fluorescent mixture can be illustrated through the sim¬ 
plified Xe-TMA diagram in Fig. 21. The relevant reaction 
rates are kq (Xe* quenching by TMA), kp (Penning trans¬ 
fer rate), (fluorescent transfer rate) and ksQ (TMA self¬ 
quenching rate), of which the first three are unknown. In 
this discussion we interpret these rate constants as effec¬ 
tive parameters, keeping in mind that more complex be¬ 
haviours may ensue, e.g., with fluorescence sequentially fol¬ 
lowing Penning. 

First, it must be noted that a lower bound to the pri¬ 
mary scintillation expected from a Xenon-TMA mixture 
can be derived from studies at low pressure, by consid¬ 
ering only direct TMA excitation/scintillation and self¬ 
quenching. Hence, we conservatively omit the presence of 
Xe*, Xe^ —»TMA*(3s) transfers (fcp ~ 0), and the quench¬ 
ing of TMA excited states by Xenon species, following [29]. 
Such a bound results in: 




Si(TMA) = 


W sc t 3s {P ■ c- k S Q + l/r 3s ) 


(13) 


as a function of concentration (c) and pressure (P). By con¬ 
sidering for TMA a typical value W sc = 30 ± 10 eV (~ 
5.5 x E 3s ), the measured lifetime of the 3s state (r 3s = 44 
ns) and the TMA self-quenching rate (ksQ = 1-08 ns" 1 
bar" 1 ), evaluation of eq. 13 leads to Si yields up to 100’s 
of ph/MeV, depending on pressure [45]. Additional scin¬ 
tillation hangs on the ability of TMA to quickly dissipate, 
collisionally, a fairly large amount of internal excess energy 
(~ 2 eV) following the Xe*, Xe?S —>TMA* transfers; relax¬ 
ation to the 3s state may in this way happen while avoiding 
the dissociative pathways [29]. Although with limited effi¬ 
ciency, this transfer process is strongly hinted by the mea¬ 
surements in [40]. 


TMA 


Xe 


Xe, 



Fig. 21. Main processes involved in a ‘Penning-Fluorescent’ 
Xenon-TMA mixture. Approximate potential surfaces are repre¬ 
sented as a function of the principal coordinate (for details the reader 
is referred to [30,46,47]). Energy levels display the correct scale: for 
TMA E 3s = 5.460 eV, E 3pxy = 6.161 eV, E 3pz = 6.229 eV, adia¬ 
batic ionization potential IP = 7.74 eV, vertical ionization potencial 
IPy = 8.44 eV, Xsdnt ~ 300 nm; for Xe E 3 p 2 = 8.31 eV, E 3Pl = 8.43 
eV; for Xe 2 E 3 ^+ = 7.92 eV, E ls + = 8.05 eV, A s dnt ~ 172 nm. 
Only the most obvious path-ways are shown (arrows): Penning and 
fluorescent transfer rates (characterized by kp and kp) as well as 
Xe* quenching rates in the presence of TMA (&q) remain unknown 
at the moment. An estimate of the Penning probability (rp) as a 
function of kp, kQ and the rate of excimer formation ( k exc ) is given 
later in text. 

To the extent that energy transfers occur with the Xe* 
states, Penning reactions are energetically viable: 8 


functional behaviour of the Penning probability rp has 
been derived in [41], as: 


rp = 


cPkp 


(1 - c) 2 P 2 k exc + cP(k P + k Q ) + 1/ t* 


(15) 


where r* represents an effective life-time of the state from 
which Penning takes place. Maximum Penning transfer 
rates in the range 10-30% are hence obtained, decreasing 
with pressure. 

7.3. Outlook 

It has been shown that the use of trimethylamine as 
VUV-quencher in Xenon results in a 20-fold reduction of 
the transverse diffusion relative to pure Xenon. Besides al¬ 
lowing the stable operation of microbulk Micromegas at 
10 bar, Xe-TMA exhibits substantial primary scintillation 
above 250 nm, at the level of ~ 100’s ph/MeV [40]. The en¬ 
ergy resolution obtained for 0.511 MeV q’s is at least x2 
better than the one obtained earlier with wire chambers, 
with Xe-TMA reducing electron diffusion in factors of ~ 2 
in all 3 space dimensions relative to, e.g., Xe-CH 4 . 

From the above considerations it may seem that, 
if present radiopurity levels of microbulk Micromegas 
(0.4 Bq/m 2 ) can be reduced to the levels expected from 
the raw materials alone, re-sizing the present concept ‘as 
is’ could greatly improve on earlier limits obtained for 
/3/30rz-decay in high pressure Xenon in [8]. In this scenario, 
the optimum pressure is undoubtedly an essential vari¬ 
able to be considered. Gains of 2000 were shown for small 
sensors at 5 bar in [22] as compared to 400 at 10 bar. A 
maximum gain of 200 could be achieved for the present 
700cm 2 -readout system, limiting the energy resolution. 

Despite the positive comparison with earlier attempts in 
charge-mode readouts, it remains unclear to this collabora¬ 
tion how such a technological approach could outperform 
an electro-luminescence option, specially if more conven¬ 
tional low-diffusion mixtures at tolerable VUV-quenching 
rates do exist, as it seems [53]. On the other hand, although 
much more challenging, the studies performed on Xe-TMA 
show that the existence of a magic ‘Penning-Fluorescent’ 
mixture with high Penning and (wavelength-shifted) fluo¬ 
rescence yields cannot be discarded out of hand. Work in 
this direction will continue. 


Xe* + TMA Xe + TMA+ + e~ 


(14) 


8 . Conclusions 


The increased charge collection associated to the presence 
of Penning transfers can be a beneficial agent, reducing 
the Fano factor [49] and increasing the ionization response 
of the admixture; in the case of charge multiplication the 
avalanche fluctuations are reduced at any given gain [50]. 
By comparing Magboltz simulations with data obtained 
from charge amplification in Micromegas, an approximate 


8 Other channels are generally present [48], but they do not change 
the general argument followed here. 


A Time Projection Chamber (NEXT-MM) housing 1.1 
kg of Xenon in its fiducial volume and instrumented with 
a microbulk Micromegas readout plane (manufactured ‘a 
la’ GEM) has been operated for 100 live-days in a 10 bar 
Xenon/trimethylamine ‘Penning-Fluorescent’ gas mixture 
(99/1). The TPC ran continuously and without the need 
of dedicated shifts, showing a Micromegas damage rate at 
the halting of operation extrapolable to 1%/year. A seem¬ 
ingly reducible 5% damage was observed during system 
























commissioning. The loss of exposure due to the instabilities 
introduced by the observed rate of (non-damaging) micro¬ 
discharges in the readout was quantified to be less than 
0.05%. 

We have detailed the calibration procedures and fully 
characterized the TPC response with Xenon X-rays. The 
energy resolution (FWHM) ultimately achieved on the 
full fiducial volume was 14.6% at 30 keV, with contri¬ 
butions from the limited S/N, sampling frequency and 
non-uniformity of the readout plane explaining the dete¬ 
rioration with respect to results obtained earlier in small 
setups (9%). At the Qpp of 136 Xe, and assuming ^fe scal¬ 
ing, the calorimetric response to 0.511 MeV and 1.275 MeV 
electron tracks anticipates FWHM energy resolutions from 
3.2%(for the best sector) to 3.9%(full TPC). We have pro¬ 
posed concrete ways to further mitigate the damage rates 
and to improve the readout. They could bring the energy 
resolution to the level obtained for X-rays (1.6% at Qpp), 
or even to the intrinsic value obtained for small sensors 
(1.0% at Qpp). 

A simplified topological analysis carried out with the sys¬ 
tem’s native voxelizationof 8nnnx8mmx 1.2mm ( xxyxz) 
indicates a strong enhancement of the mip-to-blob identifi¬ 
cation capabilities in the z-direction to the extent of com¬ 
pletely decoupling the nrip and blob energy distributions 
for a sub-sample of 600 selected tracks. Due to the very low 
electron diffusion measured for the mixture (at the scale 
of 1 mm for 1 m drift), and the easiness at increasing the 
readout granularity, the technology introduced thus offers 
the possibility of mm-accurate true-3D reconstruction of 
MeV-electron tracks on large detection volumes and at high 
pressure. In particular, it was shown that the electron end¬ 
point or ‘blob’ (2mm-cr at 10 bar) could be well resolved 
with little ambiguity. 

Taking the performance of the NEXT-MM TPC at 
face value, the scaling of the present system ‘as is’ could 
supersede earlier attempts to find the f3/30v decay with 
wire-chambers. In order to be competitive with a tech¬ 
nology based on electroluminescence, a solution to the 
‘T 0 -problem’, as well as an enhanced calorimetric response, 
seem to be needed. Also, by virtue of their low material 
budget and radiopure materials, Microbulk micromegas 
hold a promise for ultra-high radiopurity, therefore efforts 
to improve the present bounds of around 0.4 Bq/m 2 for 
the 238 U and 232 Th chains should continue. 
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